Magnetic-field induced quantum-size cascades in superconducting nanowires 
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In high-quality nanowires, quantum confinement of the transverse electron motion splits the band 
of single-electron states in a series of subbands. This changes in a qualitative way the scenario of the 
magnetic-field induced superconductor-to-normal transition. We numerically solve the Bogoliubov- 
de Gennes equations for a clean metallic cylindrical nanowire at zero temperature in a parallel 
magnetic field and find that for diameters D < 10 -f- 15 nm, this transition occurs as a cascade 
of subsequent jumps in the order parameter (this is opposed to the smooth second-order phase 
transition in the mesoscopic regime). Each jump is associated with the depairing of electrons in 
one of the single-electron subbands. As a set of subbands contribute to the order parameter, the 
depairing process occurs as a cascade of jumps. We find pronounced quantum-size oscillations of 
the critical magnetic field with giant resonant enhancements. In addition to these orbital effects, 
the paramagnetic breakdown of Cooper pairing also contributes but only for smaller diameters, i. 
e., D < 5nm. 

PACS numbers: 74.78.-w, 74.78.Na 



I. INTRODUCTION 

High-quality superconducting nanostructures as, e.g., 
single-crystal Sn nanowires 1 , polycrystalline (but made 
of strongly coupled grains) Al nanowires 2 ^ and single- 
crystalline atomically uniform Pb nanofilm£ 5 »i£i£i£ have 
recently been fabricated. It was possible to minimize 
the disorder such that the electron mean free path was 
about or larger than the specimen thickness.—^ In this 
case the scattering on nonmagnetic imperfections influ- 
ences only the electron motion parallel to the wire/film, 
while the perpendicular electron motion is governed by 
the transverse-size quantization. Indeed, photoemission 
spectra of ultrathin single-crystal Pb films showed clear 
signatures of the splitting of the electron band into a 
series of subbands due to the transverse-size quantiza- 
tion^ In the presence of minimal disorder the so-called 
Anderson theorem^ (see, also, discussion in Ref. [Io|) con- 
trols the effect of nonmagnetic impurities. Thus, one 
can expect that the study of a clean system with quan- 
tized transverse electron motion can capture important 
issues concerning the impact of quantum confinement 
on the superconducting characteristics in high-quality 
nanowires /nanofilms. 

The single-electron subbands appearing due to the 
transverse quantization, move in energy with changing 
specimen thickness. When the bottom of a subband 
passes through the Fermi surface, the density of single- 
electron states at the Fermi level increases abruptly. This 
results in size-dependent superconducting resonances 11 
and in quantum-size oscillations of the superconduct- 
ing properties as function of the thickness. Recently 
such quantum-size oscillations in the critical tempera- 
ture T c of superconducting Pb nanofilms were observed 
at a high level of experimental precision and sophistica- 
tion^ Quantum-size superconducting resonances were 



shown to be responsible for an increase of the supercon- 
ducting transition temperature in Al and Sn nanowires 
with decreasing thickness 1 ^. 

The transverse quantization of the electron motion 
should strongly influence the superconducting-to-normal 
phase transition driven by a magnetic field in such high- 
quality nanowires/nanofilms. In the present paper we 
limit ourselves to nanowires in a parallel magnetic field 
and ignore the vortex formation because vortices cannot 
nucleate in very thin superconducting wires. 

According to the Ginzburg-Landau (GL) theor y 13 : 14 , 
the critical magnetic field is expected to increase as 
1/D in the Meissner state, with D the diameter of the 
mesoscopic wire. Furthermore, the superconducting-to- 
normal phase transition in a magnetic field is of second 
order for such mesoscopic wires while being of first order 
in bulk (for type I superconductors)^. It is a general char- 
acteristic of the GL theory that this transition becomes of 
second order in mesoscopic specimens^ 1 ^. Recent calcu- 
lations based on the Bogoliubov-de Gennes (BdG) equa- 
tions 16 for wires with diameters 20 -j- 200 nm confirmed 
the GL result and revealed a smooth superconducting- 
to-normal transition in a parallel magnetic field at any 
temperature below T c . This is in agreement with recent 
experimental data on Sn 1 -! 3 - and Zni£ nanorods. Hence, 
one may conclude that effects of the transverse quantiza- 
tion of the electron motion are not significant for metallic 
superconducting wires with width larger than 20 nm. 

In the present paper we show that the situation 
changes dramatically for smaller widths. Our analysis is 
based on a numerical self-consistent solution of the BdG 
equations for a clean cylindrical metallic nanowire. We 
predict that at zero temperature the superconducting-to- 
normal transition driven by a magnetic field parallel to 
the nanowire, occurs as a cascade of jumps in the order 
parameter (with clear signatures of hysteretic behavior) 
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for diameters D < 10 -r- 15 ran. This qualitative change is 
accompanied by pronounced quantum-size oscillations of 
the critical magnetic field with large enhancements at the 
points of the superconducting resonances. In addition 
to these orbital effects, we found that Pauli paramag- 
netism can also contribute but its role is only significant 
for smaller diameters, i.e., D < 5nm. 



II. BOGOLIUBOV-DE GENNES EQUATIONS 

In the clean limit the BdG equations 9 read 

E n u n {r) =H e u n (r) + A(r)v n (r), (la) 



A*(r) u n {v)-H* e v n {v), 



(lb) 



where A(r) stands for the superconducting order param- 
eter (* for complex conjugate), E n is the quasiparticle 
energy, u n (r) and v n (r) are the particle- like and hole- like 
wave functions. The single-electron Hamiltonian appear- 
ing in Eqs. (flaj) and (jlbj) is given by 



H e = 



1 
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with m e the electron band mass (can be set to the 
free-electron mass without loss of generality), and Ep 
the Fermi level. The BdG should be solved in a self- 
consistent manner, together with the self-consistency re- 
lation 



A(r) 



9 u n( T ) v U r ){ 1 



2/n), 



(3) 



with g the coupling constant and /„ = f(E n ) the Fermi 
function 9 -. 

An important issue is the range of the states included 
in the sum in Eq. ([3]). The usual prescription concerns 
the quasiparticles with positive energies E n . At the same 
time the corresponding single-electron energy should 
be located in the Debye window, |£ n | < Hujjj with the 
Debye frequency and 



d 3 r u;(r)if e tt n (r) +v*(r)H e v n (i 



(4) 



However, in the presence of a magnetic field, this pre- 
scription is modified: H e \A=o is used rather than H e in 
Eq. It is well-known that the selection |£„| < hwn 
appears as a result of the delta-function approximation 
for the effective electron-electron interaction. Such an ap- 
proximation neglects a complex structure of the Fourier 
transform of the pair interaction. The problem is cured 
by the well-known cut-off in the canonical-momentum 
space. Such a cut-off results in the above selection rule for 
£ n with H e replaced by H c \a=o (see, for instance, Refs.H 
and [l6h. Second, the requirement of positive quasipar- 
ticle energies has to be weakened in the presence of a 



magnetic field. Namely, one needs to include the states 
having positive quasiparticle energies only at zero mag- 
netic field. This allows one to investigate also the regime 
of gapless superconductivity when the presence of quasi- 
particles with negative energies manifests the depairing 
reconstruction of the ground state (see Eq. © below and 
Appendix [K§ . 

Due to transverse quantum confinement we set 



Un(r)\ 



r£S 



M r )\ 



res 



= 



(5) 



on the wire surface. Periodic boundary conditions are 
used along the nanowire. Screening of the external mag- 
netic field can be neglected for narrow wires. Then, for a 
constant magnetic field parallel to the nanocylinder, H\\, 
it is convenient to use the well-known Coulomb gauge. 
Thus, for cylindrical wires we have A(r) = A(p) with 
p, ip, z the cylindrical coordinates (below the order pa- 
rameter is chosen as a real quantity) . The set of relevant 
quantum numbers is n = {j,m, fc}, with j the quantum 
number associated with p, m the azimuthal quantum 
number, and k the wave vector of the quasi-free electron 
motion along the nanowire. In this case the particle-like 
and hole-like wave functions can be represented as 



u n (r) = u jmk (p) 



v„(r) = v jmk (p) 



2tt x/X' 



/2¥ sfL" 



(6a) 



(6b) 



with L the length of the nanowire. Inserting Eqs. (|6aj) 
and (|6bj) into the BdG equations (jlap and (jlbj) and using 
an expansion in terms of the Bessel functions (see details 
in Ref. ll8T). the problem is reduced to the diagonalization 
of a matrix. 



III. DISCUSSION OF NUMERICAL RESULTS 

A. Resonances in the critical magnetic field 

At a superconducting resonance the main contribu- 
tion to the different superconducting quantities comes 
from the subband (or subbands) whose bottom passes 
through the Fermi surface. For cylindrical wires, the sub- 
bands with the same |m| are degenerate for H\\ — and, 
hence, any size-dependent resonant enhancement of the 
order parameter (e.g., the energy gap and the critical 
temperature) can be specified by the set (j, \m\) in the 
absence of a magnetic field. Due to quantum-size os- 
cillations in the pair-condensation energy, we get corre- 
sponding oscillations in the critical magnetic field whose 
resonant enhancements can also be labeled by (j, |m|). 
Figure [TJa) shows the critical field H c \\ calculated self- 
consistently from Eqs. (fTa[) and (| lb[) at zero temperature 
(T = 0) for an aluminum nanocylinder with diameter 
D. Note that H c u is set as the magnetic field above 
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FIG. 1: (Color online) (a) Critical parallel magnetic field versus the nanowire diameter D, and spatially averaged order 
parameter A as function of H\\ for the resonant diameters (b) D — 1.94 nm [governed by (j, \m\) = (1,0)] and (c) D = 
3.1 nm [governed by (j, \m\) = (1, 1)]. 



which the spatially averaged order parameter A drops 
below O.OlAbuik, with Ab u ik the bulk gap. Here, we con- 
sider as an example Al and take 9 hw^ = 32.31 meV and 
gN(0) = 0.18, where N(0) stands for the bulk density 
of states. For this choice Abuik = 0.25 meV. The ef- 
fective Fermi level is set to Ep — 0.9 eV, which is used 
together with the BdG equations within the parabolic 
band approximation*-. As seen, H c \\ exhibits huge en- 
hancements as compared to the bulk critical magnetic 
field -ff c , bulk = 0.01 T (to simplify our discussion, we 
show first the results for extremely narrow quantum 
wires). Resonances in H c u are found to be very de- 
pendent on D and \m\. The states with large |m| are 
more strongly influenced by H\\ and, so, the resonances 
in H c \\ governed by large |m| are, as a rule, less pro- 
nounced. In contrast, the resonances controlled by m = 
are very stable. For instance, a superconducting solu- 
tion to Eqs. ([Taj) and (|lb[) exists at D = 1.94 nm [the 
resonance associated with (j, \m\) — (1,0)] even for an 
abnormally large magnetic field of about 1000 T. Simi- 
lar behavior is found for the resonance at D = 3.21 nm 
with (j,m) = (2,0). Note that in Fig. Q]two neighbor- 
ing resonances with (j,m) = (2,0) [D — 3.21 nm] and 
(j, \m\) = (1,5) [D = 3.28 nm] merge and result in one 
profound increase in H c ,||- 



B. Quantum-size cascades 

Figures HJb) andQJc) show two typical examples (m = 
and |m| ^ 0) of how the spatially averaged order pa- 
rameter A depends on Hu. To discuss these results, we 
remark that the quasiparticle energies can be well ap- 
proximated by 





>7 



A 2 , 



(7) 



where £j m k is the single-electron energy given by Eq. (0]) 
(at -ffn = 0), Hb stands for the Bohr magneton and 



Aj- mfc = dp p A(p) \u jmk (p)\ 2 + \v jmk {p)Y 



(8) 



the averaged value of the order parameter as seen by 
jmk— quasiparticles (R = D/2). Equation §8§ can be de- 
rived within Anderson's approximate solution of the BdG 
equations 2 ^. This approximate solution implies that the 
particle-like and hole-like wave functions are chosen to 
be proportional to the eigenfunctions of H e (for details, 
see Appendix [AJ . Note that the dependence of Aj m k on 



k is found to be negligible: Aj m k = A jm (see Eq. (|A6|) ). 
As follows from Eq. ([7]), quasiparticles with in > are 
moved down in energy by Hu. Each time when a quasi- 
particle branch specified by a positive m touches zero, a 
jump in A occurs. When a branch controlling a resonant 
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FIG. 2: (Color online) The quasiparticle energies Ej m k versus k for the four relevant branches (j,m) = (0, 0), (0, ±1), (0, ±2) 
and (1,0) at (a) Hu = and (b) Hu = 33. 5T for the resonant diameter D = 1.94nm. In this case there are two nonzero 
solutions of the BdG equations for 32.9 T < H\\ < 33.5 T [see the left-side inset in Fig. D^b)], the upper one disappears at 
Hi I = 33.5 T [panel (b)] when the quasiparticle branch (0,2) touches zero, (c) The quantity mink Ej m k versus H\\ for the 
different quasiparticle branches at D = 1.94 run. 



enhancement approaches zero, A jumps down to zero and 
the superconducting solution disappears (see Fig. QJc)). 
Other quasiparticle branches arc less important (due to 
a smaller density of states) and are responsible for small 
(sometimes almost insignificant) jumps in A. In par- 
ticular, at D — 1.94 nm the first small jump in A (see 
Fig. [lib)) is located at ifi = 33.5 T. Here the branch 
with j = 0,m = 2 touches zero (see Fig. [5]). The in- 
sets in Fig. [ljb) show details of jumps in A. As seen, 
there are clear signatures of hysteretic behavior: in the 
vicinity of any jump the BdG equations has two possible 
solutions. 

To properly clarify details of the hysteretic behav- 
ior, we performed a numerical analysis for sufficiently 
large values of the unit-cell length L, controlling peri- 
odic boundary conditions in the longitudinal direction. 
In particular, the limit L — > oo can be approached only 
when L > 10 20 ^m (L/D > 10 5 ). For m_ = the 
last term in Eq. ([7]) is "switched off' and, so, A exhibits 
only a sequence of small jumps for the resonances gov- 
erned by m = [see Fig. QJb)]. For any quasiparticle 
branch an energy gap A% m) (see Fig. [2) can be intro- 
duced, and the total excitation energy gap is defined as 
Ae = minA^ . Stress that in general, Aj m ^ A^ m \ 
only at if 1 1 = we have A^ = min Aj m . Thus, a jump 

in A appears when one of A^, m ^ becomes zero. In par- 
ticular, the left-side inset in Fig. [TJb) shows that there 
exist two solutions in the interval from if| = 32.94 T 
to if 1 1 = 33.5 T (the first jump in the order parameter 
as a function of Hu). For the upper solution we have 



A^ = A0.2 > that decreases linearly with if| until 
touching zero at Hu = 33.5 T [see the quasiparticle en- 
ergies corresponding to the upper solution and given in 
Figs. (Ha) andEHJb) for if|| = and H\\ = 33.5 T, respec- 
tively]. For the lower solution A^?' 2 -* = and, so, the 

gapless regime is realized with = A^' 2 ' = 0. For 
more detail, Fig. [UJc) shows how mint Ejmk varies with 
if 1 1 for D = 1.94 nm. As seen, each relevant quasiparti- 
cle branch exhibits signatures of two small jumps [corre- 
sponding to the jumps in the order parameter given in 
Fig. [2f b)] . After the first jump min k -Eo,2,fc < (see the 
inset) and, so, A^?' 2 ' = 0. After the second jump (for 
Hu > 55.85 T) min k -Eo,i,fc becomes negative and, hence, 
we get Ae = A^' 1 "* = A^' 2 -* = 0. In the near vicin- 
ity of the second jump, for 55.7 T < H\\ < 55.85 T, we 
again find two nonzero solutions for the BdG equations: 
A^?' 1 '* ^ for the upper solution (except of the edge point 

ifn = 55.85 nm) and A^?' 1 '' = for the lower one. Above, 
we discussed only numerical results for the resonances. 
The same conclusions hold for the off-resonant points. 
However, the eventual jump to zero in A at if || = if Ct n 
is, of course, much less pronounced in this case. 

Now the question arises: what is the physics underlying 
these cascades of jumps in the order parameter? A jump 
appears when one of the relevant quasiparticle branches 
touches zero. From this point on, such a branch "sup- 
plies" the system with states having negative quasiparti- 
cle energies [see the discussion about Eq. © in Sec. ITT] . 
For such quasiparticles /„ = 1 at zero temperature or, in 
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FIG. 3: (Color online) The upper (a,b) and lower (c,d) solutions of the BdG equations at H\\ — 33.5 T (the resonant diameter 
D — 1.94 nm): (a) and (c), the order parameter A(p); (b) and (d), the corresponding contribution of the different relevant 
single-electron subbands. 



other words, these quasiparticles survive even at T = 0. 
It means that we face a reconstruction of the ground 
state. To have a feeling about such a reconstruction, let 
us consider the multiband BCS ansatz for the ground- 
state wave function (see Appendix [5] Eq. (|A9|) V This 
ansatz reads 

l*> = II Wmk - V* mk ai mM al_ m> _ kl )\0), (9) 

where aj m k ^ (a^ m ^) is the creation (annihilation) op- 
erator for electrons in the state j, m, k with the z spin 
projection "f, and Uj mk and Vj mk are given by 

Ujmk = J d 3 r if* jmk (r)u jmk (r), (10a) 

Vjmk = J d 3 r (p* jmk (r)vj mk (r), (10b) 

with (pj m k(r) being the eigenfunction of H e (the term 
oc A 2 (r) can be neglected, see, for instance, Ref. [ill ). 

/2~ a ■ e lmi P e ikz 



where J m (x) is the mth order Bessel function, and ay m 
is its j— th zero. When a quasiparticle with a negative 
energy appears at T = (say, with the quantum numbers 
j', ml ', k', |), the ground state given by Eq. ([9]) should be 
abandoned in favor of 

Ti/.m'.fc/.fl*) = 

°5', m ',*M II ( U U - V* mk al mM al_ m> _ ktl )\0), (12) 

j m' k 

where 7J, m , k , » stands for the quasiparticle creation op- 
erator, 

7j' , m ' , k ' ,1 = Uj'm'k' a j' jm \ k '1 + Vj'm'k' 0>j' ,-m' ,-k' ,|- 

As seen, Eq. (fT2|) differs from Eq. ([9]) due to the sector 
j',m',k': in Eq. (|12[) we simply have the single-electron 
creation operator rather than the Cooper-pair correlation 
term including the product a - m k ^aj _ m ^ . Therefore, 
the reconstruction mentioned above is due to depairing of 
electrons. For instance, as seen from Figs.[5{b) and[^c), 
the quasiparticle branch with j = 0, m = 2 touches zero 
at = 33. 5T and, at higher magnetic fields, acquires 
negative energies. This gives rise to the depairing of elec- 
trons in the single-electron subband j = 0, m = 2, which 
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results in the drop of the order parameter [see the left- 
side inset in Fig.QJb)]. Note that such a drop occurs not 
only due to a decay of the Cooper pairs in the subband 
j = 0,m = 2. Throughout the self-consistency relation 
([3]) , such a decay influences and reduces the contributions 
of all other subbands. However, the binding energies of 
the Cooper pairs in these subbands are somewhat re- 
duced rather than the deparing of electrons occurs. In 
Fig. [3] the order parameter is plotted together with the 
contributions of different single-electron subbands for the 
upper [(a) and (b)] and lower [(c) and (d)] solutions of 
the BdG equations at Hu = 33.5 T and D = 1.94nm. 
Comparing panels (a) and (c), we find that the order 
parameter decreases slightly by a few percent, which re- 
sults in a small jump of A in Fig. QJb) (the left-side in- 
set). From Figs. [3f b) and[3Jd), we can see that all the 
subband contributions are also reduced by a few percent 
when passing from the upper to the lower solution, except 
for j = 0, to = 2. For j = 0, m = 2 we have a significant 
drop by a factor of 1.5, which is a manifestation of elec- 
tron depairing. In a quasi-one-dimensional system there 
is a set of single-electron subbands contributing to the 
order parameter, and, so, the depairing process occurs as 
a quantum-size cascade of jumps. 
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C. Effect of thickness 

In the previous subsection, for the sake of simplicity, 
we considered extremely small diameters. So the ques- 
tion arises about the effect of thickness. In Fig. 0Ja) A 
is plotted as a function of Hu and D for larger diame- 
ters, i.e., D = 4 -4- 6nm. We see that the quantum-size 
oscillations in H c u are correlated (as to the positions 
of the resonances) with the corresponding oscillations in 
A. However, contrary to the A-resonances, amplitudes of 
resonant enhancements in H c u are mainly determined by 
1 777 1 . The most profound increases in H c u correspond to 
?77 = and appear at D = 4.55 and 5.9 nm. Signatures of 
jumps in A can again be observed (see, also, the contour 
plot given in Fig.QJb)). For instance, at D = 4.22 nm the 
averaged order parameter A jumps from a value about 
2 meV down to zero at -Ef|| = i? c ,|| = 6 T. At D = 4.77 nm 
a jump of about 1 meV occurs at Hu = H c u — 4T. For 
the off-resonant values of D we also have jumps in A 
but less pronounced. Note that the resonant enhance- 
ments in the superconducting condensate governed by 
j = 3, m = and j — 4, to = are very stable against 
iJ| | but decay significantly faster as compared to the sit- 
uation of smaller diameters. 

The number of relevant single-electron subbands scales 
as cx D 2 , which results in complex patterns of the hys- 
teretic behavior accompanying the jumps in A at larger 
diameters. An example of such a complex pattern is 
shown in Fig. [5l where details of the first jump in A(Hu) 
are given for D = 4.26 nm (the half-decay point of the 
resonance appearing at D = 4.22 nm). In this case 
there are two hysteretic loops. The larger loop is re- 



FIG. 4: (Color online) (a) Averaged order parameter A as 
a function of H\ | and D for diameters 4-^6 nm and (b) the 
contour plot of this function. The dashed curve in (b) shows 
the GL result for the critical magnetic field. 

alized for 3.02 T < H\\ < 3.16 T (see panel (a)). Sur- 
prisingly, it includes a smaller hysteretic loop arising for 
3.08 T < < 3.1 T. In this magnetic-field range there 
exist three solutions of the BdG equations. Low- lying 
quasiparticle energies for each of these solutions are given 
in Figs.E[b)-(d) for H\\ = 3.09 T. As seen, all the quasi- 
particle energies are positive for the upper solution (panel 
(b)), which is the gap regime and A^ = Ao,7 > 0. For 
the middle solution (panel (c)) we have mink-Eo,7,fc < 0, 
and, so, Aoj = 0. This is a signature of the depair- 
ing of electrons in the subband with j = 0, m = 7. For 
the lowest solution (panel (d)) the decay of the Cooper 
pairs occurs in the two single-electron subbands with the 
quantum numbers j = 0, m = 7 and j = 0, m = 6. For 
both the middle and lowest solution negative quasiparti- 
cle energies make a contribution to the problem, which 
is typical of the gapless regime. 

Note that the Ginzburg-Landau (GL) theory is not 
able to give the found quantum-size cascades and the 
quantum-size oscillations in the critical magnetic field 
(due to the absence of quantum confinement in the GL 
formalism). When using a simplified estimate based on 

the GL formula 13 ' 14 H^ L) = 8.0AiJ Cibulk /Z) (with A 
the magnetic penetration depth) together with the zero- 
temperature expectations A ~ 50 nm and H c = 0.01 T 
for Al in the clean limit^), we obtain the dashed curve 
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FIG. 5: (Color online) (a) Multi-hysteretic behavior at the first jump in A as a function of H\\ for D = 4.26 nm, and the lowest 
quasiparticle energies for three solutions of the BdG equations at H\\ = 3.09 T: (b) the upper with (b) A = 1.01 meV, (c) the 
middle with A = 0.88 meV, and (d) the lower with A = 0.83 meV. 



in Fig. BJb) which gives roughly the averaged trend for 
H c \\ found with the BdG formalism. 

For thicker mesoscopic wires with D > 20 nm, the role 
of any given quasiparticle branch becomes much less sig- 
nificant, and quantum-size oscillations in the supercon- 
ducting properties are strongly reduced. In this regime 
we recover the smooth superconducting-to-normal transi- 
tion in agreement with the previous theoretical results^ 
and recent experimental observation o 1 ' 17 ? 21 . 



D. Pauli paramagnetism 

We remark that in the current approach we neglected 
Pauli paramagnetism entirely and included only orbital 
effects. This is justified when the paramagnetic (Pauli) 
limiting field 22 

A E (H = 0) 
n P = 7x 

is larger than the orbital values of H c n (note that A ps 
As at zero magnetic field). From Fig. |4jb) one can es- 
timate that Hp k, 23 T (12, 14, 16, 11, 9 and 14 T) ver- 
sus iJ Ci || w 6T (9, 4, 11, 3, 2 and 4T) at D = 4.22 nm 
(4.45, 4.77, 5.2, 5.33, 5.68 and 5.85 nm). As seen, Pauli 



paramagnetism is only crucial for the resonances gov- 
erned by m = 0, i.e., at D = 4.55 and 5.88 nm, and it 
can produce some minor corrections to the resonances 
governed by \m\ = 1 (see, for instance, D — 5.2 nm) 
and by \m\ = 2 (see, for example, D = 4.45 nm). How- 
ever, most of the resonant enhancements for D > 5 nm 
are produced by the states with |m| > 2 (the larger the 
diameter, the smaller the relative number of resonances 
labeled by \m\ < 2). 

Thus, our numerical results are not very sensitive to 
the spin-magnetic interaction for D > 5 nm, whereas sig- 
natures of jumps in A are observed up to D « 10-^15nm. 
In particular, Fig. [6] shows A (a) and A^ (b) versus H\\ 
at D = 10.86 nm (H P = 3.9 T) and D = 11.24 nm (H P = 
4.5 T). The energy gap decays as a set of lines with differ- 
ent slopes, which reflects the linear dependence of Ej m k 
on f/)| in Eq. ([7]). It is remarkable that only jumps to zero 
in A are clearly seen in Fig. (Ha): a cascade of preceding 
small jumps has nearly collapsed into a continuous curve. 



IV. CONCLUDING REMARKS 

The quantization of the transverse electron motion 
in high-quality nanowires results in the splitting of the 
single-electron band into a series of subbands. Based on 
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FIG. 6: (Color online) (a) Spatially averaged order parameter 
A and (b) the energy gap Ae versus H\\ for the resonant 
diameters D = 10.86 nm and D — 11.24 nm. 



electron subbands. Thus, the same conclusions should 
hold for superconducting high-quality films (but not for 
nanograms where the orbital effects are known to be 
negligible, see, for instance, Refs. and l2~ih . It is 
well-know n 25 ' 26 that for ultrathin films the paramagnetic 
breakdown of the Cooper pairing results in a first-order 
superconducting-to-normal transition driven by a paral- 
lel magnetic field (provided that the effect of the spin- 
orbital scattering is not very significant and the temper- 
ature is close to zero). We expect that the quantum- 
size cascades can precede this paramagnetic breakdown. 
Fluctuations in thickness can somewhat destroy the cas- 
cades, and, so, atomically uniform high-quality nanofilms 
should be used to observe the orbital effects predicted in 
this paper. 
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APPENDIX A: ANDERSON'S APPROXIMATE 
SOLUTION 



a numerical solution of the Bogoliubov-de Gennes equa- 
tions for a clean metallic nanocylinder, we showed that 
such a splitting leads to important qualitative changes in 
the interplay of superconductivity and magnetic field in 
nanowires with diameters < lO-f-15 nm. At zero tempera- 
ture the superconducting-to-normal transition driven by 
a parallel magnetic field occurs as a cascade of jumps 
in the order parameter (a second-order phase transition 
is realized for mesoscopic wires). At the same time the 
critical magnetic field exhibits quantum-size oscillations 
with pronounced resonant enhancements. 

Our results are for nanowires with uniform cross sec- 
tion along the wire. Real samples will exhibit inevitable 
cross-section fluctuations that will smooth the quantum- 
size oscillations of superconducting properties, resulting 
in an overall enhancement with decreasing thickness [for 
H c \\ this enhancement can follow the simple estimate 
based on the GL theory, see Fig. [3fb)]. Such a monoton- 
ical increase of T c has recently been found in Al and Sn 
nanowires^ 2 -. At present, the parallel critical magnetic 
field has been measured in Sn^i and Zn¥- wires with 
diameters down to 20 nm. These nanowires were found 
to be still in the mesoscopic regime. It is expected that 
data on ff c i| for D < 20 nm will be available in the near 
future. 

Note that, on the qualitative level, our results are not 
sensitive to the specific confining geometry, the only thing 
that is of importance, is the formation of the single- 



To have an idea about the validity of Eqs. ([7]) and ©, 
it is instructive to consider Anderson's approximate so- 
lution to the BdG equations^. The main assumption is 
that Wn(r) and v n (r) are proportional to the cigenfunc- 
tion of H e given by Eq. (fPTj) (the term oc A 2 (r) in H e 
can be ignored), 



u n (r) = U„ ip n (r), v n (r) = V n y>„(r), 



(Al) 



with n = {j,m,k}. Note that Uj m k and Vj m k are the 
same as m Eqs. flTUaj) and (flUbf . Inserting Eqs. (fAT|) 
into Eqs. (fTa|) and (flb|) . we recast the BdG equations 
into 



E 



jmk U jmk 



^jmk—^B'm-H^ Ujmk+AjmkVjmk, (A2a) 



^■jmkUjmk 



£,jmk+HBmH\\ 



V jmk , (A2b) 



jmk (t ne order parameter is chosen as 



where A jm fc = A 

real) is given by Eq. ([S]), [Ib stands for the Bohr magne 
ton and £,j m k (the single-electron energy at H\\ = 0, see 
the discussion next to Eq. ||3J) is of the form 



2m P 



R 2 



(A3) 



with otj m the jth zero of the mth order Bessel function. 
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Equations (|A2a[) and (|A2b[) have a nontrivial solution 
only when the relevant determinant is zero, 



Ejmk — S,jmk + —Aj mk 

— A jink Ej m k + ijmk + HBfnH.\\ 

which leads to 



E jmk = ±^£j mk + Aj mk - p B mH\ | , (A4) 

where the + sign corresponds to the physical solution. 
This explains Eq. |(7J) used for the interpretations of our 
numerical results in Sec. IIII1 Taking into account the 
normalization condition (Uj m k and Vj m k are real) 



U 



jmk ' v jmk 



y , ., 



(A5) 



together with Eqs. ([8]) and (|A1[) , one can find that Aj m j. 
does not depend on k [see our discussion after Eq. |8])] , 



A 



jink 



W = S dppJU^p) (A6) 



Now, for a given Aj m , Eqs. (|A2a[) an d ()A2b[) can be 
solved analytically, which results in (for the physical 

Ejmk) 



II 2 = 
jmk 



v 2 , 

jmk 



£,jmk 



1 - 



UjmkVjmk 



\J (jink + 


A 2 ' 


^jmk 




If 2 + 

y ^>jmk ' 


A 2 ^ 







A 2 

jm 



(A7a) 
(A7b) 

(A7c) 



It is worth noting that the magnetic field is not present 
explicitly in Eqs. (|A7ap and (|A7b|l . and Uj rn k and Vj rn k 
depend on H\\ only through Aj m . Equations (|A1|) . (|A7a|) 
and (|A7b[) make it possible to rewrite Eq. §5§ in the form 
of the following BCS-like self-consistency equation: 

a A jm tanh{(3 E jmk /%) , 

/_\j/ m f - 2 vj'm'.jm i j V^-^-V 



jmk 



2je . + a 2 



with (3 the inverse temperature and the pair-interaction 
matrix element 



Vv 



2g 

ttR 2 L 



1 1 

JdppJ m ,(- 



'■j'm' 

R 



R 



-p)- 



The summation in Eq. (( A8|) is over the physical states 
with £j m k being in the Debye window, |^m/c| < Fhod- 

Note that Eqs. ()A1[) is exact only when A(p) = const, 
which is not the case in the presence of quantum con- 
finement. However, one can expect that Anderson's ap- 
proximation is good enough when the Cooper pairing 
of electrons from different subbands is negligible, i.e., 
for narrow wires with a strong impact of the trans- 
verse quantization. This expectation is in agreement 
with our numerical results revealing that Anderson's 
approximation is accurate within a few percent when 
D < 5 -j- 10 urn. In particular, according to Eq. (|A8[) . 
the superconducting order parameter is constant at zero 
temperature until quasiparticles with negative energies 
appear: tanh( / 9£'/2) — > 1 for (3 — ► oo when E > 0, 
whereas tanh(/3i?/2) — > —1 in the opposite case. This 
explains why A given in Figs. [TJ O and [5] is practically 
independent of H\\ before the gapless regime. 

As mentioned above, Anderson's solution is a good ap- 
proximation when the Cooper-pairing of electrons from 
different subbands plays a minor role. So, Anderson's 
prescription given by Eq. (|A1|) is equivalent to the multi- 
band BCS model whose Hamiltonian can be written as 
(n = {j,m,k}) 



n a 

+ 7J y ] ^ ] Vjm J'in'< 1 n rT a n-o- a n'-o- a n'o-i (A9) 



with h — {j, — m, — k} and a the electron spin projec- 
tion. Comparing Eq. (|A9|) with the bulk reduced BCS 
Hamiltonian, one can easily generalize the well-known 
BSC ansatz for the bulk ground-state wave function to 
the multiband ansatz given by Eq. © in Sec. IIIII 



1 M. L. Tian, J. G. Wang, J. Snyder, J. Kurtz, Y. Liu, P. 
Schiffer, T. E. Mallouk, and M. H. W. Chan, Appl. Phys. 
Lett. 83, 1620 (2003); M. Tian, J. Wang, J. S. Kurtz, Y. 
Liu, M. H. W. Chan, T. S. Mayer, and T. E. Mallouk, 
Phys. Rev. B 71, 104521 (2005). 

2 M. Zgirski, K.-P. Riikonen, V. Touboltsev, and K. Aru- 
tyunov, Nano Lett. 5, 1029 (2005). 

3 L. Jankovic, D. Gournis, P. N. Trikalitis, I. Arfaoui, T. 



Cren, P. Rudolf, M.-H. Sage, T. T. M. Palstra, B. Kooi, J. 
De Hosson, M. A. Karakassides, K. Dimos, A. Moukarika, 
and T. Bakas, Nano Lett. 6, 1131 (2006). 

4 F. Altomare, A. M. Chang, M. R. Melloch, Y. Hong, and 
C. W. Tu, Phys. Rev. Lett. 97, 017001 (2006). 

5 Y. Guo, Y.-F. Zhang, X.-Y. Bao, T.-Z. Han, Z. Tang, L.-X. 
Zhang, W.-G. Zhu, E.G. Wang, Q. Niu, Z. Q. Qiu, J.-F. 
Jia, Z.-X. Zhao, and Q. K. Xue, Science 306, 1915 (2004). 



10 



6 D. Eom, S. Qin, M.-Y. Chou, and C. K. Shih, Phys. Rev. 
Lett. 96, 027005 (2006). 

7 M. M. Ozer, J. R. Thompson, and H. H. Weitering, Nature 
Physics 2, 173 (2006). 

8 M. M. Ozer, J. R. Thompson, and H. H. Weitering, Phys. 
Rev. B 74, 235427 (2006). 

9 P. G. de Gennes, Superconductivity of Metals and Alloys 
(W. A. Benjamin, New York, 1966). 

10 M. A. Skvortsov and M. V. Feigel'man, Phys. Rev. Lett. 
95, 057002 (2005). 

11 J. M. Blatt and C. J. Thompson, Phys. Rev. Lett. 10, 332 
(1963). 

12 A. A. Shancnko, M. D. Croitoru, M. Zgirski, F. M. Peeters, 
and K. Arutyunov, Phys. Rev. B 74, 052502 (2006). 

13 V. P. Silin, Zh. Eksp. Teor. Fiz. 21, 1330 (1951) (in Rus- 
sian) . 

14 O. S. Lutes, Phys. Rev. 105, 1451 (1957). 

15 V. A. Schweigert and F. M. Peeters, Phys. Rev. B 57, 
13817 (1998). 

16 J. E. Han and V. H. Crespi, Phys. Rev. B 69, 214526 
(2004). 



J. S. Kurtz, R. R. Johnson, M. Tian, N. Kumar, Z. Ma, 
S. Xu, and M. H. W. Chan, Phys. Rev. Lett. 98, 247001 
(2007). 

A. A. Shanenko and M. D. Croitoru, Phys. Rev. B 73, 
012510 (2006). 

A. A. Shanenko, M. D. Croitoru, and F. M. Peeters, Euro- 
phys. Lett. 76, 498 (2006); ibid, Phys. Rev. B 75, 014519 
(2007). 

P. W. Anderson, J. Phys. Chem. Solids 11, 26 (1959). 
T. G. Sorop and L. J. de Jongh, Phys. Rev. B 75, 014510 
(2007). 

D. Saint- James, G. Sarma, and E. J. Thomas, Type II Su- 
perconductivity (Pergamon, New York, 1969). 
J. von Delft, Ann. Phys. (Leipzig) 10, 219 (2001). 
V. N. Gladilin, V. M. Fomin, and J. T. Devreese, Phys. 
Rev. B 70, 144506 (2004). 
K. Maki, Phys. Rev. 148, 362 (1966). 
P. M. Tedrow, R. Meservey, and B. B. Schwartz, Phys. 
Rev. Lett. 24, 1004 (1970). 



